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ABSTRACT 

We compute mass outflow rates from advective accretion disks around 
compact objects, such as neutron stars and black holes. These computations, for 
the first time, are done using combinations of exact transonic inflow and outflow 
solutions which may or may not form standing shock waves. Assuming that 
the bulk of the outflow is from the effective boundary layers of these objects, 
we find that the ratio of the outflow rate and inflow rate varies anywhere from 
a few percent to even close to a hundred percent (i.e., close to disk evacuation 
case) depending on the initial parameters of the disk, the degree of compression 
of matter near the centrifugal barrier, and the polytropic index of the flow. 
Our result, in general, matches with the outflow rates obtained through a fully 
time-dependent numerical simulation. In some region of the parameter space 
when the standing shock does not form, our results indicate that the disk may 
be evacuated and may produce quiescence states. 



Subject headings: black holes — neutron stars- stars: accretion — - stars:winds 
— shocks waves — jets — AGN — quasars 

Submitted to ApJ in April 1998; Still trying to teach the referee the meaning of centrifugal 
barrier and the funnel wall. 



1. Introduction 

One of the signatures of activities around compact objects is the presence of jets 
and outflows. Outflows carry away angular momentum from the accretion disk and are 
partially responsible for the accretion itself. Active galaxies and quasars are supposed to 
harbor black holes at their centers and at the same time produce cosmic radio jets through 
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which immense amount of matter and energy are ejected out of the core of the galaxies 
(See, Begelman, Blandford & Rees, 1984 and Chakrabarti, 1996a for reviews). Similarly, 
micro-quasars have also been discovered very recently where outflows are formed from 
stellar black hole candidates (Mirabel & Rodriguez, 1994). Many of these outflows show 
superluminal motions which are probably due to magnetic effects. With hydrodynamic 
effects alone, which we are employing in this paper, one should not be able to accelerate 
the flow more than the initial sound velocity. A well known stellar object SS433 (which is 
believed to be a neutron star), where pure hydrodynamic effects may be operating, produces 
outflows with a speed roughly one-third of the speed of light. 

There are several models in the literature which study the origin, formation and 
collimation of these outflows. Difference between stellar outflows and outflows from these 
systems is that the outflows in these systems have to form out of the inflowing material only. 
This is because black holes and neutron stars have no atmospheres of their own. The models 
present in the literature are roughly of three types. The first type of solutions confine 
themselves to the jet properties only, completely decoupled from the internal properties 
of accretion disks. They study the effects of hydrodynamic or magneto-hydrodynamic 
pressures on the origin of jets (e.g., Blandford & Payne 1982; Ferrari et al., 1985; Fukue 
1987; Contopoulus 1985; Chakrabarti, 1990). In the second type, efforts are made to 
correlate the internal disk structure with that of the outflow using both hydrodynamic (e.g., 
Chakrabarti 1986) and magnetohydrodynamic considerations (Konigl 1989; Chakrabarti & 
Bhaskaran 1992). In the third type, numerical simulations are carried out to actually see 
how matter is deflected from the equatorial plane towards the axis (e.g., Hawley, Smarr 
& Wilson, 1984; Eggum, Katz & Coroniti, 1985; Molteni, Lanzafame & Chakrabarti, 
1994; Molteni, Ryu & Chakrabarti, 1996; Ryu, Chakrabarti & Molteni, 1997). From the 
analytical front, although the wind type solutions and accretion type solutions come out 
of the same set of governing equations (Chakrabarti 1990), there were no attempt to find 
connections among them. As a result, the estimation of the outflow rate from the inflow 
rate has always been impossible. On the other hand, the mass outflow rate of the normal 
stars are calculated very accurately from the stellar luminosity. Theory of radiatively driven 
winds seems to be very well understood (e.g., Castor, Abott & Klein, 1975). The simplicity 
of black holes and neutron stars lie in the fact that they do not have atmospheres. But the 
disks surrounding them have, and similar method as employed in stellar atmospheres should 
be applicable to the disks. Our approach in this paper is precisely this. We first determine 
the properties of the rotating inflow and outflow and identify solutions to connect them. In 
this manner we self-consistently determine the mass outflow rates. 

Before we present our results, we describe basic properties of the rotating inflow and 
outflow. A rotating inflow with a specific angular momentum \{x) entering into a black 
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hole will have angular momentum A ~ constant close to the black hole for any moderate 
viscous stress. This is because the viscous time scale is generally much longer compared to 
the infall time scale and even though at the outer edges the angular momentum distribution 
may be Keplerian or even super-Keplerian matter would be highly sub-Keplerian close to 
the black hole. This is because the flow has to enter through the horizon with velocity 
of light and presence of this large inertial (ram) force, in addition to usual gravitational 
and centrifugal forces, makes the flow sub-Keplerian. Almost constant angular momentum 
produces a very strong centrifugal force \ 2 /x 3 which increases much faster compared to 
the gravitational force ~ GM/x 2 and becomes comparable at around x ~ l 2 /GM, or, 
r cb ~ 2 A 2 where r and A are x and /, written in units of R g = 2GM/c 2 and R g c = 2GM/c 
respectively. The subscript cb under r stands for the centrifugal barrier. (In the rest of the 
paper, we use R g as the length unit, c is the unit of velocity, and the mass of the black 
hole M to the unit of mass.) Here, (actually, a little farther out, due to thermal pressure) 
matter starts piling up and produces the centrifugal pressure supported boundary layer 
(CENBOL). Further close to the black hole, the gravity always wins and matter enters the 
horizon supersonically after passing through a sonic point. CENBOL may or may not have 
a sharp boundary, depending on whether standing shocks form or not. Generally speaking, 
in a polytropic flow, if the polytropic index 7 > 1.5, then shocks do not form and if 7 < 1.5, 
only a region of the parameter space forms the shock (C96b). In any case, the CENBOL 
forms. In this region the flow becomes hotter and denser and for all practical purposes 
behaves as the stellar atmosphere so far as the formation of outflows are concerned. Inflows 
on neutron stars behave similarly, except that the 'hard-surface' inner boundary condition 
dictates that the flow remain subsonic between the CENBOL and the surface rather than 
becoming supersonic as in the case of a black hole. In case where the shock does not form, 
regions around pressure maximum achieved just outside the inner sonic point would also 
drive the flow outwards. In the back of our mind, we have the picture of the outflow as 
obtained by Hawley, Smarr & Wilson (1984), Chakrabarti (1986) and Molteni, Ryu and 
Chakrabarti (1996), namely, that the outflow is thermally and centrifugally accelerated but 
confined by external pressure of the ambient medium. 

There are two surfaces of utmost importance in flows with angular momentum. One 
is the 'funnel wall' where the effective potential (sum of gravitational potential and the 
specific rotational energy) vanishes. In the case of a purely rotating flow, this is the 'zero 
pressure' surface. Flows cannot enter inside the funnel wall because the pressure would be 
negative. The other surface is called the 'centrifugal barrier'. This is the surface where 
the radial pressure gradient of a purely rotating flow vanishes and is located outside the 
funnel wall simply because the flow pressure is higher than zero on this surface. Flow with 
inertial pressure easily crosses this 'barrier' and either enters into a black hole or flows out 
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as winds depending on its initial parameters (detail classification of the parameter space 
is in Chakrabarti, 1989) In numerical simulations (e.g., Hawley, Smarr & Wilson, 1984; 
Molteni, Ryu and Chakrabarti, 1996) it is observed that the outflow generally hugs the 
'funnel wall' and goes out in between these two surfaces. In this paper we assume precisely 
this. However, we believe that the the cross section area could be weakly influenced by the 
pressence of the small outward radial motion of the flow. We have ignored such a correction. 

Outflow rates from accretion disks around black hole and neutron stars must be related 
to the properties of CENBOL which in turn, depend on the inflow parameters. Subsonic 
outflows originating from CENBOL would pass through sonic points and reach far distances 
as in wind solutions. Assuming free-falling conical polytropic inflow and isothermal outflows 
(as in stellar winds), it is easy to estimate the ratio of outflowing and inflowing rates 
(Chakrabarti, 1997a; 1998): 

M ™t _ D _ Qout-R -(/„-g) f 3/2 

• -UrH- e Jo 

lvl m ^~>in ^ 

where, Q out and Oj n are the solid angles of the outflow and inflow respectively, and 

_ (2n + l)R 
k ~ 2n ■ 

Here, R is the compression ratio of inflowing matter at the CENBOL and n — 1/(7 — 1) 
is the polytropic constant. When Q out ~ Q in , ~ 0.052 and 0.266 for 7 = 4/3 and 5/3 
respectively. Assuming a thin inflow and outflow of 10° conical angle, the ratio R^ becomes 
0.0045 and 0.023 respectively. 

The aim of the present paper is to compute the mass loss rate more realistically 
than what has been attempted so far. We calculate this rate as a function of the inflow 
parameters, such as specific energy and angular momentum, accretion rate, polytropic 
index etc. We explore both the polytropic and the isothermal outflows. Our conclusions 
show that the outflow rate is sensitive to the specific energy and accretion rate of the 
inflow. Specifically, when the outflow is not isothermal, outflow rate generally increases 
with the specific energy and the polytropic index 7 of the outflow, generally decreases 
with the polytropic index 7 of the inflow, but somewhat insensitive to the specific angular 
momentum A. In the case of isothermal outflow, however, mass loss rate is sensitive to the 
inflow rate, since the inflow rate decides the proton temperature of the advective region of 
the disk which in turn fixes the outflow temperature. In this case the outflow is at least 
partially temperature driven. The outflow rate is also found to be anti-correlated with 
specific angular momentum A of the flow. 

The plan of this paper is the following: In the next Section, we describe our model and 
present the governing equations for the inflow and outflow. In §3, we present the solution 
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procedure of the equations. In §4, we present results of our computations. Finally, in §5, 
we draw our conclusions. 



2. Model Description and Governing Equations 

2.1. Inflow Model 

For the sake of computation of the inflow quantities, we assume that the inflow is 
axisymmetric and thin: h(r) << r, so that the transverse velocity component could be 
ignored compared to the radial and azimuthal velocity components. We consider polytropic 
inflows in vertical equilibrium (otherwise known as 1.5 dimensional flows, see, Chakrabarti, 
1989; hereafter referred to as C89). We ignore the self-gravity of the flow and viscosity 
is assumed to be significant only at the shock, if present. We do the calculations using 
Paczyriski-Wiita (1980) potential which mimics surroundings of the Schwarzschild black 
hole. The equations (in dimensionless units) governing the inflow are: 
(a) Conservation of specific energy is given by, 

,-> n A^ I . . 

£ = - + na < + 2?- - W^T) ■ (1) 

where, u e and a e are the radial and polytropic sound velocities respectively. a e = {^jp e / Pe) 1 ^ 2 , 
p e and p e are the pressure and density of the flow. For a polytropic flow, p e = KpJ, where 
K is a constant and is a measure of entropy of the flow. Here, A is the specific angular 
momentum and n is the polytropic constant of the inflow n = (7 — I 7 is the polytropic 
index. The subscript e refers that the quantities are measured on the equatorial plane. 

Mass conservation equation, apart from a geometric constant, is given by, 

M in = u e p e rh e (r), (2) 

where h e (r) is the half-thickness of the flow at radial co-ordinate r having the following 
expression 

(3a) 

Another useful way of writing the mass inflow is to introduce an entropy dependent quantity 
M oc ^y n K n M which can be expressed as 

M = u e aM{r-l)s^ (36) 
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Where, M is really the entropy accretion rate (C89). When the shock is not present, M 
remains constant in a polytropic flow. When the shock is present, Ai will increase at the 
shock due to increase of entropy, a — (7 + l)/(7 — 1) = 2n + 1. If the centrifugal pressure 
supported shock is present, the usual Rankine-Hugoniot conditions, namely, conservations 
of mass, energy and momentum fluxes across the shock are to be taken into account (C89) 
in determining the shock locations. In presence of mass loss one must incorporate this effect 
in the shock condition (see, eq. 15 below). 

2.2. Outflow Models 

We consider two types of outflows. In ordinary stellar mass loss computations (e.g. 
Tarafdar, 1988, and references therein), the outflow is assumed to be isothermal till the 
sonic point. This assumption is probably justified, since copious photons from the stellar 
atmosphere deposit momenta on the slowly outgoing and expanding outflow and possibly 
make the flow close to isothermal. This need not be the case for outflows from compact 
sources. Centrifugal pressure supported boundary layers close to the black hole are very hot 
(close to the virial temperature) and most of the photons emitted may be swallowed by the 
black holes themselves instead of coming out of the region and depositing momentum onto 
the outflow. Thus, the outflows could be cooler than isothermal flows. In our first model, 
we choose polytropic outflows with same energy as the inflow (i.e., no energy dissipation 
between the inflow and outflow) but with a different polytropic index 7 < 7. Nevertheless, 
it may be advisable to study the isothermal outflow to find out the behavior of the extreme 
case. Thus an isothermal outflow is chosen in our second model. In each case, of course, we 
include the possibility that the inflow may or may not have standing shocks. 

It is to be noted that we chose the outflows to be propagating through the funnel wall 
even though the inflow is assumed to be thin. This may not be completely self-inconsistent. 
On the one hand, our assumption of thin inflow is for the sake of computation of the 
thermodynamic quantities only, but the flow itself need not be physically thin. Secondly, 
the funnel wall and the centrifugal barrier are purely geometric surfaces, and they exist 
anyway and the outflow could be supported even by ambient medium which may not 
necessarily be a part of the disk itself. So, we believe that our assumptions are not entirely 
unjustified. 
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2.2.1. Polytropic Outflow 
In this case, the energy conservation equation takes the form: 

« ^ / 2 A 2 1 

f = T + nV + 2^w"2(F3i) (4) 

and the mass conservation in the outflow takes the form: 

M ut = P$A(r). (5) 

Here, n! = — l)^ 1 is the polytropic constant of the outflow. The difference between eq. 
(4) and eq. (1) is that, presently, the rotational energy term contains 

, . U(r) + R(r) . . 

r m (r) = V ' 2 , (6a) 

as the mean axial distance of the flow. The expression of 3?(r), the local radius of the 
centrifugal barrier comes from balancing the centrifugal force with the gravity (Molteni, 
Ryu & Chakrabarti, 1996), i.e., 

A 2 £(r) 



K 3 (r) 2r(r - If 

We thus obtain, 

r 2X 2 r(r - l) s 



(66) 



(7a) 



And the expression for R(r), the local radius of the funnel wall, comes from vanishing of 
total effective potential, i.e., 

1 A 2 

<W(r) = -^— yy + ^y = ° 

i?(r) = A[(r-l)] 1 / 2 (76) 

The difference between eq. (5) and eq. (2) is that the area functions are different. Here, 
A(r) is the area between the centrifugal barrier and the funnel wall (see introduction for the 
motivation). This is computed with the assumption that the outflow is external pressure 
supported, i.e., the centrifugal barrier is in pressure balanced with the ambient medium. 
Matter, if pushed hard enough, can cross centrifugal barrier in black hole accretion (the 
reason why rapidly rotating matter can enter into a black hole in the first place). An 
outward thermal force (such as provided by the CENBOL) in between the funnel wall and 
the centrifugal barrier causes the flow to come out. Thus the cross section of the outflow is, 



A{r) = n[^ 2 (r)-R 2 (r)}. 



(8) 
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The outflow angular momentum A is chosen to be the same as in the inflow, i.e., no viscous 
dissipation is assumed to be present in the inner region of the flow close to a black hole. 
Considering that viscous time scales are longer compared to the inflow time scale, it may be 
a good assumption in the disk, but it may not be a very good assumption for the outflows 
which are slow prior to the acceleration and are therefore, prone to viscous transport of 
angular momentum. Such detailed study has not been attempted here particularly because 
we know very little about the viscous processes taking place in the pre-jet flow. Therefore, 
we concentrate only those cases where the specific angular momentum is roughly constant 
when inflowing matter becomes part of the outflow, although some estimates of the change 
in R,h is provided when the average angular momentum of the outflow is lower. Detailed 
study of the outflow rates in presence of viscosity and magnetic field is in progress and 
would be presented elsewhere. 



2.2.2. Isothermal Outflow 

The integration of the radial momentum equation yields an equation similar to the 
energy equation (eq. 4): 

■&■ 2 X 2 1 

+ C 2 lnp + ? = Constant (9) 

2 2r m (r) 2 2(r - 1) 1 ' 

In this case the thermal energy term is different, behaving logarithmically. The constant 
sound speed of the outflow is C s . The mass conservation equation remains the same: 

M out = pS iso A{r). (10) 

Here, the area function remains the same above. A subscript iso of velocity d is kept 
to distinguish from the velocity in polytropic case. This is to indicate the velocities are 
measured here using completely different assumptions. 

In both the models of the outflow, we assume that the flow is primarily radial. Thus 
the ^-component of the velocity is ignored (fig «•&). 



3. Procedure to solve for disks and outflows simultaneously 

Before we go into the details, a general understanding of the transonic flows around 
a black hole is essential. In Chakrabarti (1989, 1996c), all the solution topologies of 
the polytropic flow in pseudo- Newtonian geometry has been provided (e.g., see, Fig. 3 



- 9- 



Chakrabarti, 1996c) In regions I and O of the parameter space the flow has only one sonic 
point. Matter with positive energy at a large distance must pass through that point before 
entering into the black hole supersonically. In regions SA and SW shocks may form in 
accretion and winds respectively, but no shocks are expected in winds and accretions if 
parameters are chosen from these branches. In NSW and NSA, two saddle type sonic 
points exist, but no steady shock solutions are possible. 

Suppose that matter first enters through the outer sonic point and passes through a 
shock. At the shock, part of the incoming matter, having higher entropy density is likely 
to return back as winds through a sonic point, other than the one it just entered. Thus 
a combination of topologies, one from the region SA and the other from the region O is 
required to obtain a full solution. In the absence of the shocks, the flow is likely to bounce 
back at the pressure maximum of the inflow and since the outflow would be heated by 
photons, and thus have a smaller polytropic constant, the flow would leave the system 
through an outer sonic point different from that of the incoming solution. Thus finding a 
complete self-consistent solution boils down to finding the outer sonic point of the outflow 
and the mass flux through it. Below we present the list of parameters used in both of our 
models and briefly describe the procedure to obtain a satisfactory solution. 



3.1. Polytropic Outflow 

We assume that 

(a) In this case, a very little amount of total energy is assumed to be lost in each bundle 
of matter as it leaves the disk and joins the jet. The specific energy S remains fixed 
throughout the flow trajectory as it moves from the disk to the jet. 

(b) Very little viscosity is present in the flow except at the place where the shock forms, so 
that the specific angular momentum A is constant in both inflows and outflows close to the 
black hole. At the shock, entropy is generated and hence the outflow is of higher entropy 
for the same specific energy. 

(c) The polytropic index of the inflow (7) and outflow (7 Q ) are free parameters and in 
general, 7 D < 7, because of heating effect of the outflow (e.g., due to the momentum 
deposition coming out of the disk surface). In reality 7 D is directly related to the heating 
and cooling processes of the outflow. When M; n is high, heating of outflow by photon 
momentum deposition is higher, and therefore 7 D — > 1. 

Thus a supply of parameters S, A, 7 and 7 D make a self-consistent computation of 
possible when the shock is present. When the shock is absent, the compression ratio of the 
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gas at the pressure maximum between the inflow and outflow R comp is supplied as a free 
parameter, since it may be otherwise very difficult to compute satisfactorily. In the presence 
of shocks, such problems do not arise as the compression ratio is obtained self-consistently. 

The following procedure is adopted to obtain a complete solution: 



(a) From eqs. (1) and (2) we derive an expression for the derivative, 

A 2 na 2 5r — 3 1 



du _ r 3 a r(r — 1) 2(r - l) 2 , , 

~dr ~ 2ru? ' ^ ' 

u 



era 



At the sonic point, the numerator and denominator separately vanish, and give rise to the 
so-called sonic point conditions: 

1 A 2 



a c = 2( V 1)2 n T(? (12a) 
a{r c - l)r c 

n(5r c — 3) 

[2n . 
u c = \ — a c (12b) 
V a 

where, the subscript c represents the quantities at the sonic point. The derivative of 
the flow at the sonic point is computed using the L'Hospital's rule. Using fourth order 
Runge-Kutta method i?(r) and a(r) are computed along the flow till the position where the 
Rankine-Hugoniot condition is satisfied (if shocks form) and from there on the sub-sonic 
branch is integrated for the accretion as usual (C89). With the known j , £ and A, one can 
compute the location of the outflow sonic point from eqs. (4) and (5), 

a 2 dA[r) A 2 dr m (r) 1 



d$ _ A 2 (r) dr r m 3 (r) dr 2(r - l) 2 , . 

from where the sonic point conditions at the outflow sonic point r co obtained are given by, 
a 2 G dA(r) A 2 dr m (r) 1 



A 2 co (r) dr r rnco 3 (r) dr 2(r co - l) 2 

and 



(14a) 



$co = a co . (146) 
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At the outer sonic point, the derivative of d is computed using the L'Hospital's rule and the 
Runge-Kutta method is used to integrate towards the black hole to compute the velocity of 
the outflow at the shock location. The density of the outflow at the shock is computed by 
distributing the post-shock dense matter of the disk into spherical shell of An solid angle. 
The outflow rate is then computed using eq. (5). 

It is to be noted that when the outflows are produced, one cannot use the usual 
Rankine-Hugoniot relations at the shock location, since mass flux is no longer conserved in 
accretion, but part of it is lost in the outflow. Accordingly, we use, 

M + = (l-R lh )M. (15) 

where, the subscripts + and — denote the pre- and post-shock values respectively. Since due 
to the loss of matter in the post-shock region, the post-shock pressure goes down, the shock 
recedes backward for the same value of incoming energy, angular momentum & polytropic 
index. The combination of three changes, namely, the increase in the cross-sectional area 
of the outflow and the launching velocity of the outflow and the decrease in the post-shock 
density decides whether the net outflow rate would increased or decreased than from the 
case when the exact Rankine-Hugoniot relation was used. 

In the case where the shocks do not form, the procedure is a bit different. It is 
assumed that the maximum amount of matter comes out from the place of the disk where 
the thermal pressure of the inflow attains its maximum. The expression for the polytropic 
pressure for the inflow in vertical equilibrium is, 

2(n+l)/iV. 

y ' 7 (!+™).M 

This is maximized and the outflow is assumed to have the same quasi-conical shape with 
annular cross-section Air) between the funnel wall and the centrifugal barrier as already 
defined. In the absence of shocks the compression ratio of the flow between the incoming 
flow and outgoing flow at the pressure maximum cannot be computed self-consistently 
unlike the case when the shock was present. Thus this ratio is chosen freely. We take the 
guidance for this number from what was obtained in the case when shocks are formed. 
However, in this case even when the mass loss takes place, the location of the pressure 
maximum remains unchanged. Since the compression ratio R comp is a free parameter, 
remains unchanged for a given R comv - Let us assume that /i_ is the actual mass inflow 
rate and it is same before and after the pressure maximum had the mass loss rate been 
negligible. Let /i + be the mass inflow rate after the pressure maximum, when the loss due 
to outflow is taken into account. Then, by definition, /i_ = M out + /i + and R^ = M out /fi_. 
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Thus the actual ratio of the mass outflow rate and the mass inflow rate, when the mass loss 
is taken into consideration is given by, 

Mout Rrh (in \ 

~7T~ = 1 — r~- (17a) 

However, this static consideration is valid only when < 1. Otherwise, we must have, 

_ dM disk 

i.e., 



+ fl- = fi + + M out 



" /< (R,„ 1) • (176) 



dt 

Here, M rfisfc is the instantaneous mass of the disk. Since R^ > 1, the disk has to evacuate. 
These cases hint that the assumptions of the steady solution breaks down completely and 
the solutions may become become highly time dependent. 



3.2. Isothermal Outflow 

We assume that 

(a) The outflow has exactly the same temperature as that of the post-shock flow, but the 
energy is not conserved as matter goes from disk to the jet. In other words the outflow is 
kept in a thermal bath of temperature as that of the post-shock flow. 

(b) Same as (b) of §3.1. 

(c) The post-shock proton temperature is determined from the inflow accretion rate M in 
using the consideration of Comptonization of the advective region. The procedure to 
compute typical proton temperature as a function of the incoming accretion rate has been 
adopted from Chakrabarti (1997b). 

(d) The polytropic index of the inflow can be varied but that of the outflow is always unity. 

Thus a supply of parameters £, A and 7 makes a self-consistent computation of Rm 
possible when the shock is present. When the shock is absent, the compression ratio of the 
gas at the pressure maximum between the inflow and the outflow R comp is supplied as a free 
parameter exactly as in the polytropic case. 

The following procedure is adopted to obtain a complete solution: 
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(a) From eqs. (9) and (10) we derive an expression for the derivative, 

C 2 dA(r) A 2 dr m (r) 1 



M\ _ Mr) dr r m 3(r) dr 2(r c - l) 2 
dr lis °~ „ C 2 • 1 J 
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At the sonic point, the numerator and denominator separately vanish, and give rise to the 
so-called sonic point conditions: 

Cs 2 dA(r) X 2 dr m {r) 1 

A co (r) dr lc ° + r m J(r) dr |c ° 2(r co — l) 2 ~~ ^ 

and 

Ac = C„ (196) 

where, the subscript co represents the quantities at the sonic point of the outflow. The 
derivative of the flow at the sonic point is computed using the L'Hospital's rule. The 
procedure is otherwise similar to those mentioned in the polytropic case and we do not 
repeat them here. 



4. Results 

4.1. Polytropic outflow coming from the post-shock accretion disk 

Figure 1 shows a typical solution which combines the accretion and the outflow. The 
input parameters are S = 0.0005, A = 1.75 and 7 = 4/3 corresponding to relativistic 
inflow. The solid curve with an arrow represents the pre-shock region of the inflow and 
the long-dashed curve represents the post-shock inflow which enters the black hole after 
passing through the inner sonic point (I). The solid vertical line at X s3 (the leftmost 
vertical transition) with double arrow represents the shock transition obtained with exact 
Rankine-Hugoniot condition (i.e., with no mass loss). The actual shock location obtained 
with modified Rankine-Hugoniot condition (eq. 15) is farther out from the original location 
X s3 . Three vertical lines connected with the corresponding dotted curves represent three 
outflow solutions for the parameters 7 Q = 1.3 (top), 1.15 (middle) and 1.05 (bottom). The 
outflow branches shown pass through the corresponding sonic points. It is evident from the 
figure that the outflow moves along solution curves completely different from that of the 
'wind solution' of the inflow which passes through the outer sonic point 'O'. The mass loss 
ratio in these cases are 0.256, 0.159 and 0.085 respectively. Figure 2 shows the ratio 
Rjh as 7 is varied. Only the range of 7 and energy for which the shock- solution is present 
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is shown here. The general conclusion is that as 7 Q is increased the ratio is also increased 
non-linearly. When the inflow rate is very low, due to paucity of the photons, the outflow is 
not heated very much and 7 D remains higher. The reverse is true when the accretion rate 
is higher. Thus, effectively, the ratio R m is going up with the decrease in M{ n . In passing 
we remark that with the variation in the inflow angular momentum, A, the result does not 
change significantly, and R m changes only by a couple of percentage at the most. 

In Fig. 3a, we show the variation of the ratio R m of the mass outflow rate and inflow 
rate as a function of the shock-strength (dotted) M_/M + (Here, M_ and M + are the 
Mach numbers of the pre- and post-shock flows respectively), the compression ratio (solid) 

(Here, £_ and S + are the vertically integrated matter densities in the pre- and 
post- shock flows respectively), and the stable shock location (dashed) X s3 (in the notation 
of C89). Other parameters are A = 1.75 and 7 D = 1.05. Note that the ratio R m does not 
peak near the strongest shocks! Shocks are stronger when they are located closer to the 
black hole, i.e., for smaller energies. It is interesting to note that the general shape of R m 
variation roughly agrees with analytical results of Chakrabarti (1997a; 1998) where it was 
shown that for a finite outflow compression ratio must be non-zero and R m peaks at an 
intermediate compression ratio. The non-monotonic behavior is more clearly seen in lowest 
curve of Fig. 3b where R m is plotted as a function of the specific energy £ (along x-axis) 
and 7 (marked on each curve). Specific angular momentum is chosen to be A = 1.75 
as before. The tendency of the peak in R m is primarily because as £ is increased, the 
shock location is increased which generally increases the outflowing area A(r) at the shock 
location. However, the density of the outflow at the shock, as well as the velocity of the 
outflow at the shock increases. The outflow rate, which is a product of these quantities, 
thus shows a peak. For the sake of comparison, we present the results for 7 D = 1.05 (dashed 
curve) when the Rankine-Hugoniot relation was not corrected by eq. (15). The result 
generally remains the same because of two competing effects: decrease in post-shock density 
and increase in the area from the the outflow is launched (i.e., area between the black hole 
and the shock) as well as the launching velocity of the jet at the shock. 

To have a better insight of the behavior of the outflow we plot in Fig. 4 R m as a 
function of the polytropic index of the incoming flow (7) for 7 D = 1.1, £ = 0.002 and 
A = 1.75. The range of 7 shown is the range for which shock forms in the flow. We also plot 
the variation of injection velocity "&mj, injection density pi n j and area A(r) of the outflow 
at the location where the outflow leaves the disk. The incoming accretion rate has been 
chosen to be 0.3 (in units of the Eddington rate). These quantities are scaled from the 
corresponding dimensionless quantities as ~&i n j — > 0.W in j, p in j — > 10 22 p in j and A — > 10 _4 ^4 
respectively in order to bring them in the same scale. With the increase in 7, the shock 
location is increased, and therefore the cross-sectional area of the outflow goes up. The 
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injection velocity goes up (albeit very slowly) as the shock recedes, since the injection 
surface (CENBOL) comes closer to the outflow sonic point. However, the density goes down 
(gas is less denser). This anti-correlation is reflected in the peak of R^- 

So far, we assumed that the specific angular momentum of the outflow is exactly the 
same as that of the inflow, while in reality it could be different due to presence of viscosity. 
In the outflow, a major source of viscosity is the radiative viscosity whose coefficient is, 

4 « T4 2-1 

i] = — cm sec (20) 

This could be significant, since the temperature of the outflow is high, but the density is 
low. Assuming that the angular momentum distribution reaches a steady state inside the 
jet (Chakrabarti, 1986 and references therein), 

lj = CjR n > (21) 

where Cj and rij are constants, the vanishing condition of the azimuthal velocity on the axis 
requires that rij > 1 inside the jet. The matter distribution in the rotationally dominant 
region of the 'pre-jet' is computed by integrating Euler equation. It is easy to show that the 
'hollow' jet thus produced carry most of the matter and angular momentum in the outer 
layers of the jet (Chakrabarti, 1986). In other words, the average angular momentum of 
the outflow away from the base may remain roughly constant even in presence of viscosity. 
This is to be contrasted with the disk, where matter is more dense towards the centre while 
more angular momentum is concentrated towards the outer edge. If, however, the average 
angular momentum at the base of the outflow goes does down due to losses to ambient 
medium, by, say, a factor of two, we find that the mass loss rate is also reduced by around 
the same factor. This shows that the outflow is at least partially centrifugally driven. 

An important point to note: the ratio between the 'specific entropy measure' of the 
outflow to that of the post-shock inflow is obtained from the definitions of entropy accretion 
rate M (see, eq. 3b): 

K _ MZi 1 (I - R^- 1 ^( 7 - 7o) 7. (22) 



K + Ml' 1 V RrH ) + 7 o 

As Rm — > 1, -f^- — > 0. Thus, we expect that for a polytropic flow with shocks, a hundred 
percent outflow is impossible since the outgoing entropy must be higher. In isothermal 
outflows such simple consideration do not apply. 

If we introduce an extra radiation pressure term (with a term like T/r 2 in the radial 
force equation, where T is the contribution due to radiative process), particularly important 
for neutron stars, the outcome is significant. In the inflow, outward radiation pressure 
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weakens gravity and thus the shock is located farther out. The temperature is cooler and 
therefore the outflow rate is lower. If the term is introduced only in the outflow, the effect 
is not significant. 

4.2. Polytropic outflow coming from the region of the maximum pressure 

In this case, the inflow parameters are chosen from region I (see C96c) so that the 
shocks do not form. Here, the inflow passes through the inner sonic point only. The 
outflow is assumed to be originated from the regions where the polytropic inflow has a 
maximum pressure. This assumption is justified, since it is expected that winds would get 
the maximum kick at this region. Figure 5a shows a typical solution. The arrowed solid 
curve shows the inflow and the dotted arrowed curves show the outflows for 7 = 1.3 (top), 
1.1 (middle) and 1.01 (bottom). The ratio R m in these cases is given by 0.66, 0.30 and 0.09 
respectively. The specific energy and angular momentum are chosen to be S = 0.00584 and 
A = 1.8145 respectively. The pressure maximum occurs outside the inner sonic point at r p 
when the flow is still subsonic. Figure 5b shows the variation of thermal pressure of the flow 
with radial distance. The peak is clearly visible. Since the pressure maximum occurs very 
close to the black hole as compared to the location of the shock, the area of the outflow is 
smaller, but the radial velocity as well as the density of matter at the base of the outflow 
are much higher. As a result the outflow rate is exorbitantly higher compared to the shock 
case. Figure 6 shows the ratio R m as a function of 7 D for various choices of the compression 
ratio Rcomp of the outflowing gas at the pressure maximum: R comp = 2 for the rightmost 
curve and 7 for the leftmost curve. We have purposely removed the solutions with R m > 1, 
because the solution should be inherently time-dependent (see, eq. 17b) in these cases and 
a steady state approach is not supposed to be trusted completely. This is different from 
the results of §4.1, where shocks are considered, since Rm is non-monotonic in that case. It 
is also found that in some range of parameters, the very high massflow could take place 
even for smaller compression ratios especially when the sonic point of the outflow r Q is right 
outside the pressure maximum. These cases can cause runaway instabilities by rapidly 
evacuating the disk. They may be responsible for quiescent states in X-ray novae systems 
and also in some systems with massive black holes (e.g., our own galactic centre?), although 
it is not clear if such a phase of profuse mass loss has been directly observed. 

The location of maximum pressure being close to the black hole, it may be generally 
very difficult to generate the outflow from this region. Thus, it is expected that the ratio 
R m would be larger when the maximum pressure is located farther out. This is exactly 
what we see in Fig. 7, where we plot R m against the location of the pressure maximum 
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(solid curve). Secondly, if our guess that the outflow rate could be related to the pressure is 
correct, then the rate should increase as the pressure at the maximum rises. That is also 
observed in Fig. 7. We plot here R m as a function of the actual pressure at the pressure 
maximum (dotted curve). The mass loss is found to be a strongly correlated with the 
thermal pressure. Here we have multiplied non-dimensional thermal pressure by 1.5 x 10 24 
in order to bring it in the same scale. 



4.3. Isothermal outflow coming from the post-shock accretion disk 

In this case, the outflow is assumed to be isothermal. The temperature of the outflow 
is obtained from the proton temperature of the advective region of the disk. The proton 
temperature is obtained using the Comptonization, bremsstrahlung, inverse bremsstrahlung 
and Coulomb processes (Chakrabarti, 1997b and references therein). Figure 8 shows the 
effective proton temperature and the electron temperature of the post-shock advective 
region as a function of the accretion rate (in units of Eddington rate, in logarithmic scale) of 
the Keplerian component of the disk. The diagram is drawn for a black hole of mass 1OM . 
The soft X-ray luminosity for stellar mass black holes or the UV luminosity of massive black 
holes is basically dictated by the Keplerian rate of the disk. It is clear that as the accretion 
rate of the Keplerian disk is increased, the advective region gets cooler as is expected. 

In Fig. 9a, we show the ratio R m as a function of the accretion rate (in units of 
Eddington rate) of the incoming flow for a range of the specific angular momentum. In the 
low luminosity objects the ratio is larger. Angular momentum is varied from A = 1.7 (top 
curve), 1.725 (middle curve) and 1.75 (bottom curve). The specific energy is £ = 0.003. 
Here we have used the modified Rankine-Hugoniot relation as before (eq. 15). The ratio 
Rm is clearly very sensitive to the angular momentum since it changes the shock location 
rapidly and therefore changes the post-shock temperature very much. We also plot the 
outflux of angular momentum F(X) = \m in R m which has a maximum at intermediate 
accretion rates. In dimensional units, these quantities represent significant fractions of 
angular momentum of the entire disk and therefore the rotating outflow can help accretion 
processes. Curves are drawn for different A as above. In Fig. 9b, we plot the variation 
of the ratio directly with the proton temperature of the advecting region. The outflow is 
clearly thermally driven. Hotter flow produces more winds as is expected. The angular 
momentum associated with each curve is same as before. 
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4.4. Isothermal outflow coming from the region of the maximum pressure 

This case produces very similar result as in the above case, except that like Section 4.2 
the outflow rate becomes very close to a hundred percent of the inflow rate when the proton 
temperature is very high. Thus, when the accretion rate of the Keplerian flow is very small, 
the outflow rate becomes very high, close to evacuating the disk. As noted before, this may 
also be related to the quiescent state of the X-ray novae. 

5. Conclusions 

In this paper, we have computed the mass outflow rate from the advective accretion 
disks around galactic and extra-galactic black holes. Since the general physics of advective 
flows are similar around a neutron star, we believe that the conclusions may remain roughly 
similar provided the shock at X s3 forms, although the boundary layer (which corresponds 
to X s i in C89 notation; also see C96b) of the neutron star, where half of the binding energy 
could be released, may be more luminous than that of a black hole and may thus affect the 
outflow rate. We have chosen a limited number of free parameters just sufficient to describe 
the inflow and only one extra parameter for the outflow. We find that the outflow rates can 
vary from a very few percentage of the inflow rate, to as much as the inflow rate (causing 
almost complete evacuation of the accretion disk) depending on the inflow parameters. For 
the first time, it became possible to use the exact transonic solutions for both the disks and 
the winds and combine them to form a self-consistent disk-outflow system. Although we 
present results when centrifugally supported boundary layers are considered around a black 
hole, it is evident that the result is general, namely, if such a barrier is produced by other 
means, such as pre-heating (Chang & Ostriker, 1985) or by pair-plasma pressure (Kazanas 
& Ellison, 1986), outflows would also be produced (Das, in preparation). 

The basic conclusions of this paper are the followings: 

a) It is possible that most of the outflows are coming from the centrifugally supported 
boundary layer (CENBOL) of the accretion disks. 

b) The outflow rate generally increases with the proton temperature of CENBOL. In other 
words, winds are, at least partially, thermally driven. This is reflected more strongly when 
the outflow is isothermal. 

c) Even though specific angular momentum of the flow increases the size of the CENBOL, 
and one would have expected a higher mass flux in the wind, we find that the rate of 
the outflow is actually anti-correlated with the A of the inflow. On the other hand, if the 
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angular momentum of the outflow is reduced by hand, we find that the rate of the outflow 
is correlated with A of the outflow. This suggests that the outflow is partially centrifugally 
driven as well. 

d) The ratio R m is generally anti-correlated with the inflow accretion rate. That is, disks of 
lower luminosity would produce higher R m . 

e) Generally speaking, supersonic region of the inflow do not have pressure maxima. Thus, 
outflows emerge from the subsonic region of the inflow, whether the shock actually forms or 
not. 

In this paper, we assumed that the magnetic field is absent. Centrifugally driven 
winds (e.g., Blandford & Payne, 1982), electromagnetically accelerated winds (Lovelace, 
1976; Lovelace, Wang & Sulkanen, 1987; Contopoulos & Lovelace, 1994 etc.) and winds 
from Blandford & Znajek (1977) process have not been considered here. These effects 
have so far been considered in the context of a Keplerian disk and not in the context of 
sub-Keplerian flows on which we concentrate here. Secondly, whereas the entire Keplerian 
disk was assumed to participate in wind formation, here we suggest that CENBOL is the 
major source of outflowing matter. It is not unreasonable to assume that CENBOL would 
still form when magnetic fields are present (Chakrabarti, 1990) and since the Alfven speed 
is, by definition, higher compared to the sound speed, the acceleration, and therefore the 
mass outflow would also be higher than what we computed here. Such works would be 
carried out in future. 

In the literature, not many results are present which deal with exact computations of 
the mass outflow rate. Molteni, Lanzafame & Chakrabarti (1994), in their SPH simulations, 
found that the ratio could be as high as 15 — 20 per cent when the flow is steady. In 
Ryu, Chakrabarti & Molteni (1997), 10 — 15 per cent of the steady outflow is seen and 
occasionally, even 150 of the inflow is found to be ejected in non-stationary cases. Our 
result shows that high outflow rate is also possible, especially for absence of shocks and 
low luminosities. In Eggum, Coroniti & Katz (1985), radiation dominated flows showed 
Rm ~ 0.004, which also agrees with our results when we consider high accretion rates (see, 
e.g., Fig. 9a). Observationally, it is very difficult to obtain the outflow rate from a real 
system, as it depends on too many uncertainties, such as filling factors and projection effects. 
In any case, with a general knowledge of the outflow rate, we can now proceed to estimate 
several important quantities. For example, it had been argued that the composition of the 
disk changes due to nucleosynthesis in accretion disks around black holes and this modified 
isotopes are deposited in the surroundings by outflows from the disks (Hogan & Applegate 
1987; Mukhopadhyay & Chakrabarti, 1998 and references therein). Similarly, it is argued 
that outflows deposit magnetic flux tubes from accretion disks into the surroundings (Daly 
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& Loeb, 1990). Thus a knowledge of outflows are essential in understanding varied physical 
phenomena in galactic environments. 

A number of possible improvements could be made on this basic work. For instance, 
the effect of radiation pressure on both the inflows and outflows are to be taken into account 
self-consistently, particularly when the neutron star winds are considered. A preliminary 
investigation with a T/r 2 force term (whose effect is to weaken gravity) suggests that for 
non-zero T in the inflow, the mass-loss rate changes significantly. This is because the shock 
location increases when T is increased. This in turn reduces the mass loss rate. On the 
other hand, the when T is non-zero in the outflow, the effect is not very high, since the 
outflow rate is generally driven by thermal effect of the disk and not the wind. Similarly, we 
see a significant reduction of the outflow when the average specific angular momentum of 
the outflow is reduced. This is expected since the outflow is partially centrifugally driven. 
This effect is stronger when the outflow is isothermal. 

An interesting situation arises when the polytropic index of the outflow is large and the 
compression ratio of the flow is also very high. In this case, the flow virtually bounces back 
as the winds and the outflow rate can be equal to the inflow rate or even higher, thereby 
evacuating the disk. In this range of parameters, most, if not all, of our assumptions may 
breakdown completely because the situation could become inherently time-dependent. It is 
possible that some of the black hole systems, including that in our own galactic centre, may 
have undergone such evacuation phase in the past and gone into quiescent phase. 

So far, we made the computations around a Schwarzschild black hole. In the case 
of a Kerr black hole (C96c), the shock locations will come closer and the outflow rates 
should become higher. Similarly magnetic field will change the sonic point locations 
significantly (C90). The mass outflow rates in these conditions are being studied and the 
results would be reported elsewhere (Das, 1998b). We made a few assumptions, some of 
which may be questionable. Our assumption of the wind formation from the post-shock 
region is very similar to the assumptions involved in computing mass outflow rate from 
stellar atmospheres. This is because the post-shock region behaves like a boundary layer. 
However, when the shocks are not formed, we assumed that the outflow comes out of the 
region of maximum pressure. Whereas, this may be a plausible assumption, it is not clear 
if strong outflows actually form from that close to the black hole. Similarly, whereas the 
assumption of isothermality in the stellar system may be justified, it is questionable if 
this can be rigorously valid in the present system since the proton temperatures are too 
high, and to maintain that high temperature till the sonic point, the disk has to deposit 
enormous energy into the winds. Nevertheless, we believe that our calculation is sufficiently 
illustrative and gives a direction which can be followed in the future for further refinements. 
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Figure Captions 

Fig. 1: Few typical solutions which combine accretion and outflow. Input parameters are 
S = 0.0005, A = 1.75 and 7 = 4/3. Solid curve with an incoming arrow represents the 
pre-shock region of the inflow and the dashed curve with an incoming arrow represents 
post-shock inflow which enters the black hole after passing through the inner sonic point 
(I). Dotted curves are the outflows for various 7 C (marked). Open circles are sonic points of 
the outflowing winds and the crossing point 'O' is the outer sonic point of the inflow. The 
leftmost shock transition {X s3 ) is obtained from unmodified Rankine-Hugoniot condition, 
while the other transitions are obtained when the mass-outflow is taken into account. 

Fig. 2: Ratio R m as outflowing polytropic index 7 is varied. Only the range of 7 for which 
the shock-solution is present is shown here. As 7 G is increased, the ratio is also increased. 
Since 7 G is generally anti-correlated with m in , this implies that R m is correlated with rrii n . 

Fig. 3: Variation of R rn as a function of (a) shock-strength M_/M + (dotted) compression 
ratio £+/£_ (solid) and the shock location X s3 (dashed) for j = 1.05 and (b) specific 
energy S for various 7 D (marked). A = 1.75 throughout. 

Fig. 4: Rm as a function of the polytropic index 7 of the inflow. The range of 7 shown is 
the range for which shock forms in the flow. Suitably scaled density, velocity and area of 
the flow (at the base of the outflow) on the disk surface are also shown. See text for details. 
Non-monotonicity in R m can be understood by the fact that the shock location, i.e., the 
area A(r in j) and velocity i?j n j of the outflow at the outflow origin go up with 7, but the 
density pi n j goes down. 

Fig. 5a: Few typical solutions for outflows forming out of an advective disk which does 
not include a standing shock wave. Incoming arrowed solid curve shows the inflow and 
the dashed arrowed curves with outgoing arrows show the outflows for 7 D = 1.3 (top), 1.1 
(middle) and 1.01 (bottom). At r p , thermal pressure of the inflow is maximum. 

Fig. 5b: Variation of thermal pressure P e of the incoming flow with radial distance. In 
a shock-free hydrodynamic flow, winds may form from the region around the pressure 
maximum. 

Fig. 6: Rm as a function of outflowing polytropic index 7 Q for various choices of the 
compression ratio R CO mp of the outflowing gas at the pressure maximum. From bottom to 
top curve, R comv = 2, 3, 4, 5, 6 & 7 respectively. 

Fig. 7: Variation of Rm (solid) as a function of the location r p of the maximum pressure 
and the non-dimensional pressure (dotted) P e (j' P ) (multiplied by 1.5 x 10 24 to bring to scale) 
are plotted. 



-25 - 



Fig. 8: Variation of effective proton (solid) and electron (dotted) temperatures in the 
advective regions of the accretion disk around a fO solar mass black hole as functions of 
the accretion rates of the Keplerian flow. As Keplerian rate increases, both protons and 
elections cool down. 

Fig. 9(a-b): (a) Variation of as functions of the Eddington rate of the Keplerian 
component of the incoming flow for a range of the specific angular momentum. In the low 
luminosity objects the ratio is larger, (b) Rm as a function of the proton temperature T p 
of the post-shock region. In (a), A = 1.7 (top curve), A = 1.725 (middle) and 1.75 (bottom 
curve). In (a) the angular momentum flux F(X) of the outflow is also shown (dashed curve). 
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